Accepted for publication in the Astrophysical Journal 

Preprint typeset using I^'T^]X style emulatcapj v. 5/2/11 



CHARACTERIZING THE COOL KOIs. V. 



KOI-256: A MUTUALLY ECLIPSING POST-COMMON ENVELOPE 
BINARY 



m 

o 

(N 

u , 

< 



6 



N. 



Philip S. Muirhead\ Andrew Vanderburg^'", Avi Shporer \ Juliette Becker'^''^'^''^^, Jonathan J. Swift\ James P 
Lloyd^, Jim Fuller^, Ming Zhao'', Sasha Hinkley'-'^, J. Sebastian Pineda\ Michael Bottom^ Andrew W. 

Howard^, Kaspar von Braun^''', Tabetha S. Boyajian*, Nicholas Law'-*, Christoph Baranec^ Reed Riddle^ A 
Ramaprakash'", Shriharsh p. Tendulkar\ Khanh Bui\ Mahesh Burse'", Pravin Chordia'", Hillol Das'", 
Richard Dekany', Sujit Punnadi'" and John Asher Johnson'-*' 

Accepted for publication in the Astrophysical Journal 

ABSTRACT 

We report that Kepler Object of Interest 256 (KOI-256) is a mutually eclipsing post-common en- 
velope binary (ePCEB), consisting of a cool white dwarf (M* = 0.592 ± O.O89Af0, i?* = 0.01345 ± 
0.00091 Rq, Teff = 7100±700 K) and an active M3 dwarf {M^ = 0.51±0.16Mo, R^, = 0.540±0.014i?Q, 
Teff = 3450 ± 50 if) with an orbital period of 1.37865 ± 0.00001 days. KOL256 is listed as hosting 
a transiting planet-candidate by Borucki et al. and Batalha et al.; here we report that the planet- 
candidate transit signal is in fact the occultation of a white dwarf as it passes behind the M dwarf. 
We combine publicly-available long- and short-cadence Kepler light curves with ground-based mea- 
surements to robustly determine the system parameters. The occultation events are readily apparent 
in the Kepler light curve, as is spin-orbit synchronization of the M dwarf, and we detect the transit 
of the white dwarf in front of the M dwarf halfway between the occultation events. The size of the 
white dwarf with respect to the Einstein ring during transit (i?Ein = 0.00473 ± 0.00055 i?©) causes the 
transit depth to be shallower than expected from pure geometry due to gravitational lensing. KOI-256 
is an old, long-period ePCEB and serves as a benchmark object for studying the evolution of binary 
star systems as well as white dwarfs themselves, thanks largely to the availability of near-continuous, 
ultra-precise Kepler photometry. 
Subject headings: binaries: eclipsing — binaries: spectroscopic — stars: abundances — stars: activity 

— stars: fundamental parameters — stars: individual (KOT256) — stars: late-type 

— stars: low-mass — stars: rotation — white dwarfs 
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1. INTRODUCTION 

Post-common-envelope binaries (PCEBs) serve as im- 
portant laboratories for studying the long-term evolu- 
tion of binary star systems (e.g. Rebassa-Mansergas et al. 
2007; Schreiber et al. 2010; Zorotovic et al. 2010, 2011; 
Morgan et al. 2012). Consisting of a co-orbiting white 
dwarf and main-sequence star, PCEBs are presumed to 
have evolved from a stage when the main-sequence star 
orbited within the envelope of the white dwarf progeni- 
tor, an asymptotic giant branch (AGB) star (e.g. Paczyn- 
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ski 1976; Willems & Kolb 2004). In this scenario, the 
main-sequence star loses angular momentum as it orbits 
within the envelope, eventually leaving a main-sequence 
star and white dwarf in a short-period orbit (e.g. Nord- 
haus & Spiegel 2012; Spiegel 2012). A PCEB will con- 
tinue to lose angular momentum via magnetic braking 
and gravitational waves and will eventually begin mass- 
transfer, resulting in a cataclysmic variable. For this rea- 
son, PCEBs are sometimes referred to as pre-cataclysmic 
variables (pre-CVs, e.g. Schreiber & Gansicke 2003). To 
date, thousands of PCEBs have been discovered spectro- 
scopically (e.g. Rebassa-Mansergas et al. 2010). 

Eclipsing PCEBs (ePCEBs) are much rarer, with sev- 
eral dozen detected to date (Drake et al. 2010; Law et al. 
2012; Parsons et al. 2012). Such systems allow for pre- 
cise measurements of the orbiting stars' physical parame- 
ters. By combining precision photometric measurements 
of the white dwarf occultation events with radial veloc- 
ity measurements, the masses, radii and temperatures of 
the component stars can be determined to relatively high 
accuracy and precision with few assumptions. 

In this paper, we show that Kepler Object of Interest 
256 (KOI-256) ^^ is an ePCEB consisting of an active 
M3 dwarf (dMe) and a cool white dwarf. KOI-256 was 
listed as hosting a transiting planet candidate in Borucki 
et al. (2011) and Batalha et al. (2013), referred to as 
KOL256.01, with an orbital period of 1.3786789 days 

1* KIC 11548140, 2MASS J 19004443-1-4933553, o=285.° 18513, 
5=-|-49.°56537 



and a fractional transit depth of 2.42 %. Batalha et al. 
(2013) note that the long-cadence planet-candidate tran- 
sit light-curve appears to be V-shaped: an indication of 
a false-positive. In the following sections, we show that 
the planet-candidate transit is in fact the occultation of 
a white dwarf as it travels behind an M dwarf, and that 
KOI-256.01 is therefore a false positive planet candidate. 
This is confirmed by large sinusoidal radial velocity vari- 
ations in the M dwarf spectrum synchronous with the oc- 
cultation ephemeris, a trapezoidal occultation light curve 
with sharp ingress and egress, and a weak transit signal 
occurring halfway between the occultation events. The 
depths of the occultation and transit are consistent with 
an eclipsing white dwarf, including a shallower transit 
depth than expected from geometry alone due to gravi- 
tation lensing. With an orbital period of greater than 1 
day, KOI-256 is a relatively long period cPCEB (Drake 
et al. 2010; Nebot Gomez-Moran et al. 2011; Law et al. 
2012; Parsons et al. 2012). As we show in Section 3.2, the 
white dwarf is also cool {Tcs ~ 7100 ± 700 K), consistent 
with an old white dwarf. 

In a previous letter, Muirhead et al. (2012) character- 
ized KOI-256 as an M3 dwarf with Toff = 3450 ± 50 
K and [M/H] = +0.31 ± 0.10 using the near-infrared 
spectroscopic indices of Rojas-Ayala et al. (2012). They 
interpolated those values onto the Dartmouth evolution- 
ary isochrones (Dotter et al. 2008; Feiden et al. 2011) 
to estimate the mass (0.43 ± O.OGM©) and radius (i?* = 
0.42 ± O.O6i?0) of the M dwarf. However, that analysis 
assumed KOI-256 was a field M dwarf, and the derived 
mass and radius did not consider the effects of PCEB 
evolution, rapid rotation or stellar activity. For exam- 
ple, Boyajian et al. (2012) found that eclipsing binary 
M dwarfs obey markedly different mass- Toff and radius- 
Teff relationships than single field M dwarfs. Stassun 
et al. (2012) investigated the effect of chromospheric ac- 
tivity on radius- Tgff relationships, due to Tcff suppression 
by the presence of spots. Rojas-Ayala et al. (2012) found 
that two nearby M dwarfs with white dwarf companions 
are significantly metal-rich compared to the solar neigh- 
borhood, leading them to speculate that the M dwarfs 
were polluted. These studies imply that the Dartmouth 
evolutionary isochrones must be used with caution when 
determining the mass and radius of a binary or active M 
dwarf based solely on its Toff and metallicity. 

In this paper we report accurate stellar parameters for 
the M dwarf in KOI-256, as well as the co-orbiting white 
dwarf, measured empirically by combining the ultra- 
precise Kepler light curve with ground-based spectro- 
scopic and adaptive optics observations. In Section 2, we 
describe the observations and the measured parameters 
used in this analysis, including a revised orbital period 
and occultation and transit ephemerides. In Section 3 
we combine the measurements to determine the physi- 
cal parameters of the system. In Section 4 we discuss 
the importance of KOI-256 for future studies of PCEBs, 
white dwarfs and exoplanets. 

2. OBSERVATIONS 

2.1. Visible-light Adaptive Optics Imaging 

We observed KOI-256 with the Robo-AO laser adap- 
tive optics and imaging system (Baranec et al. 2012) 
on the Palomar Observatory 60-inch Telescope on UT 



2012 15 July to look for contaminating sources in the 
Kepler apertures used to contract the light curves. We 
used a long-pass filter with a 600nm cut-on to more 
closely approximate the Kepler bandpass while maintain- 
ing diffraction-limited resolution. The observation con- 
sisted of a sequence of full-frame-transfer detector read- 
outs at the maximum rate of ^ 8.6 Hz for a total of 120 
s of integration time. The individual images were then 
combined using post-facto shift-and-add processing using 
KOI-256 as the tip-tilt star with 100% frame selection. 
The wavelength response of the Robo-AO visible camera 
is similar to the Kepler response, and the image accu- 
rately represents the degree of contamination by other 
stars in the Kepler apertures. We achieved high adaptive 
optics correction at visible- wavelengths: KOI-256 has a 
full-width-at-half- maximum of 0'.'12. 

Figure 1 shows the final image of KOI-256 and Figure 2 
shows a contrast curve indicating the 5cr lower-limits on 
the relative magnitudes of nearby background or contam- 
inating stars. The contrast curve was generated through 
a point spread function (PSF) subtraction algorithm cus- 
tomized for the Palomar Robo-AO imaging system, and 
based on the Locally Optimized Combination of Images 
algorithm (Lafreniere et al. 2007). Rather than using 
PSF reference images obtained sequentially in time on 
the science target, the algorithm uses reference images 
derived from the numerous other PSFs in the field of 
view. We see no evidence of contamination in the Kepler 
apertures for objects that are less than 5.0 magnitudes 
fainter than KOI256 at 1 arcsec of separation. Therefore, 
we proceed assuming all of the flux in the Kepler light 
curve is from the M dwarf-white dwarf system. 

2.2. Kepler Light Curve 

The Kepler Mission photometrically monitors targets 
in two modes: long-cadence mode and short-cadence 
mode (Koch et al. 2010). Long-cadence mode involves 
30-minute integrations and short-cadence mode involves 
1 minute integrations, both taken nearly continuously 
in the custom Kepler bandpass spanning roughly 5000 
to 9000 A. KOI-256 has publicly available long-cadence 
mode data from quarters through 11 (2009 May 
through 2012 January), and short-cadence mode data 
for quarters 6 and 7 (2010 June through 2010 Decem- 
ber) . Both datasets have several gaps in coverage due to 
quarterly rotation of the spacecraft, data downlink gaps, 
periods with the spacecraft in safe mode, and intermit- 
tent corrections to the spacecraft pointing (Batalha et al. 
2013). 

2.2.1. Evidence for Starspots and Synchronous Rotation 

We downloaded the Kepler long-cadence data from the 
Mikulski Archive for Space Telescopes (MAST) for quar- 
ters 1, 2, 3, 4, 5, 6, 7, 9, 10, and 11. No data is available 
for quarters 8 and 12 since during those quarters KOI- 
256 was located on CCD module 3, which failed per- 
manently about 20 days into quarter 4 (Batalha et al. 
2013). To clean the data from systematic trends we used 
the PyKE software (Still & Barclay 2012) and the co- 
trending vectors (Barclay et al. 2012, Section 2.3), fol- 
lowed by long-term polynomial detrending. 

Figure 3 shows a Lomb-Scargle periodogram (Lomb 
1976; Scargle 1982) of the long-cadence data while ig- 



noring in-occultation and in-transit data. The peri- 
odogram's strongest peak is at the same period as the 
orbit. Figure 3 inset shows a zoomed-in view of the 
strongest peak, which is in fact spht into several peaks. 
Using the Kepler long-cadence time stamps we estimated 
the width of a periodogram peak of a fixed period signal 
spanning the entire Kepler data time span, and found it 
is similar to the width of the individual sub-peaks. This 
suggests a slight variation in the photometric period dur- 
ing the course of Kepler monitoring. This is supported by 
analysis of different subsets of the Kepler data, marked 
by solid colored lines, where there are small differences in 
the position of the strongest peaks. The vertical dashed 
black line marks the 1.38 day orbital period, at the center 
of the group of peaks. 

Periodic, sinusoidal-like variability at the orbital pe- 
riod can arise from (1) a combination of effects induced 
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Fig. 1. — Visible-light adaptive optics image of KOI-256 by 
the Robo-AO system on the Palomar Observatory 60-inch Tele- 
scope. North is up and East is left, and the blue line rep- 
resents 10 arcsec. We overlay the Kepler apertures {red boxes) 
used to collect short-cadence data. KOI-256 happened to be 
near the edge of the Robo-AO field-of-view during integrations, 
resulting in the sharp cut-off on the left of the image. We detect 
no contaminating stars within the Kepler apertures. 
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Fig. 2. — Contrast curve showing Scr lower limits on the mag- 
nitude difference of a potential background or companion ob- 
ject for a given separation, calculated from the Robo-AO image 
above. Based on these data, we proceed assuming all of the flux 
in the Kepler light curves come from the M dwarf- white dwarf 
system. 



by the orbital motion, including beaming, ellipsoidal dis- 
tortion, reflected light (e.g., Zucker et al. 2007; Shporer 
et al. 2010; Faigler & Mazeh 2011; Shporer et al. 2011), 
or, (2) stellar chromospheric activity of a spin-orbit syn- 
chronized star. In the second scenario an active star is 
rotating at the same period as the orbit, as a result of 
tidal locking, and the non-uniform longitudinal distribu- 
tion of spots across the stellar surface results in observed 
sinusoidal-like variability at the orbital period. 

In order to investigate the origin of the photometric 
variability we visually examined the light curve. Figure 
4 shows three long-cadence light curves each spanning 
four days, and separated from each other by a few 100 
days. The occultations are seen as brief drops in relative 
flux, spanning 2-3 long-cadence measurements. A visual 
comparison of the three panels shows an evolution of the 
light curve shape. The three light curves show different 
amplitude, and a different relative phase between the oc- 
cultation and the out-of-eclipse modulation (especially 
between panels a and c). In addition, panel b shows a 
double peak modulation, at half the period seen in the 
other two panels. This evolution in the light curve shape 
is clear evidence that it originates from spots, since the 
orbital effects mentioned above should not show evolu- 
tion at the observed timescales. Moreover, given our es- 
timates of the system parameters (see Section 3), these 
orbital effects are expected to show variability of ~ 0.1%, 
an order of magnitude smaller than the observed ampli- 
tude. Finally, the splitting of the periodogram peak in 
Figure 3 has been observed for other active stars, and 
could be a sign of differential rotation (Chaplin et al. 
2013). 

Since the white dwarf brightness is only ^ 2% of the 
system total brightness in the Kepler band, it can not be 
the origin of the observed photometric variability. There- 
fore, we conclude that the dM component is spotted and 
tidally locked to the white dwarf. In Section 3 we use 
the fact that the M dwarf is spin-orbit synchronized in 
combination with the measured V sin i of the M dwarf 
and the occultation period and duration to empirically 
constrain the radius of the M dwarf. 

We note that the high signal-to-noise ratio (S/N) of the 
spot-induced variability, and the relatively short rotation 
period compared to the long time-span of the data, al- 
lows for detailed analysis of the spots pattern and their 
evolution. This can possibly lead to investigation of dif- 
ferential rotation, magnetic activity cycles, and magnetic 
interaction between the two stars. 

We also note that beyond the eclipses and quasi- 
sinusoidal variability, the light curve shows also many 
flares, with amplitudes that vary from <1% to ^10%, 
which occur with a frequency that is comparable to the 
orbital frequency. 

2.2.2. Precise Orbital Period 

We downloaded short-cadence Kepler light curve data 
of KOI-256 from MAST for quarters 6 and 7, spanning 
roughly 6 months nearly continuously. First, we mea- 
sured the orbital period of KOI-256 using the Box-fitting 
Least Squares algorithm (BLS, Kovacs et al. 2002) ap- 
plied to the short cadence data in an iterative scheme. 
For the first iteration, the parameters of the BLS al- 
gorithm were allowed to vary considerably around the 
nominal transit parameters from Batalha et al. (2013). 
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Fig. 3. — Lomb-Scargle median-normalized periodogram of KOI- 
256 long-cadence data, while ignoring in-occultation and in-transit 
data. The strongest peak is at 1.38 days, the same period as the 
orbit, and there are smaller peaks at its harmonics. This shows 
the light curve has a clear, high signal-to-noise variability with 
the same period as the orbit. The inset shows a zoom-in on the 
strongest peak (black solid line), showing it is split into a few lines. 
The orbital period is marked by a dashed black line. The colored 
lines in the inset show the periodogram for subsets of Kepler data 
(see legend). 



Successive iterations of BLS were run over smaller re- 
gions of frequency space, sampling the signal residue in 
finer detail and placing tighter constraints on the tran- 
sit duration. Our final iteration sampled the frequency 
range of 0.7247 to 0.7260 day^^ sixty times using 500 
bins across the phase folded data with the transit du- 
ration constrained to be within 0.031 and 0.032 times 
the period. The peak of the signal residue is estimated 
with two methods: fitting a Gaussian model and using a 
cubic spline interpolation between BLS samples. These 
methods return values of 1.3786570 and 1.3786526 days, 
respectively, in agreement at the level of 4 x 10~^, and 
shorter than the (Batalha et al. 2013) value by about 2 s. 
We take the mean of the two measurements as the final 
value, and we estimate the uncertainty to be 1 x 10~^ 
days, or roughly 0.8 s. Figure 5 shows the signal residue 
from the BLS analysis for a coarse frequency resolution, 
with the final high resolution analysis inlaid. The Gaus- 
sian fit and the cubic spline interpolation are shown as 
solid and dotted lines, respectively, with the final value 
of 1.3786548 days designated with the dashed line. 

2.2.3. Occultation Light Curve 

With an accurate period for the occultation events, we 
phase-folded the short-cadence data to fit an occultation 
model to the light curve. We used short-cadence data for 
the occultation fit, and neglected the long-cadence data, 
since high-cadence is required to resolve the ingress and 
egress. First, we detrended the light curve to remove the 
strong variations due to the aforementioned rotational 
modulation of the M dwarf and occasional flares. For 
each occultation event, we focused on a 3.36 hr wide 
window roughly centered around the occultation. We 
applied a running median-filter to the out-of-occultation 
data points using a box size of 14 minutes. We fitted 
a second-order polynomial to the median-filtered, out- 
of-occultation data, then divided all of the original data 
points in the 3.36 hr window by the polynomial fit, in- 



cluding the data taken during occultation. 

We then phase-folded the de-trended data using the 
revised period and binned the data into 15 s windows. 
For each bin we removed 3cr outliers and calculated the 
mean of the data points within the bin. Figure 6 shows 
the resulting de-trended and phase-folded short-cadence 
occultation light curve. 

We fitted a custom occultation model, which unlike 
planet transit models (e.g. Mandel & Agol 2002) leaves 
the ingress/egress times independent of the occultation 
depth. We assumed that the radius of the white dwarf 
was much smaller than that of the M dwarf, and mod- 
eled the shape of the occultation as a circle with a lin- 
ear limb- darkening coefficient passing behind a flat edge, 
with the occultation duration and ingress/egress dura- 
tions independent parameters. We convolved the mod- 
eled light curve with a boxcar window of length 58.8 
s to simulate the effects of the 58.8 s Kepler short ca- 
dence integrations. We fitted the model to the binned 
and phase folded data using a Levenberg-Marquardt al- 
gorithm. Varying the linear-limb darkening coefficient 
for the white dwarf had little effect on the resulting 
depth, duration or ingress/egress time relative to their 
estimated uncertainties; therefore, we chose to assume 
no limb-darkening in the white dwarf. 

The fitted parameters and corresponding \a uncertain- 
ties are summarized in Table 3. The fitted fractional oc- 
cuhation depth is 0.0238743 ± 0.000047, similar to the 
value listed for the planet-candidate signal in Batalha 
et al. (2013), but with higher precision. However, the 
fitted occultation duration is is 60.093 ± 0.037 minutes, 
and ingress and egress duration is 1.535 ± 0.067 min- 
utes. If the occultation were in fact a transiting planet, 
it would have to be Jupiter-sized in order to cause such a 
large fractional change. A transiting Jupiter-sized planet 
would result in a much larger ingress-to-occultation du- 
ration ratio: of order 1/10. The measured ingress-to- 
occultation duration ratio of ~l/40 is instead consistent 
with an Earth-sized emitting object passing behind the 
M dwarf. In the following sections we support this in- 
terpretation with a detection of the transit of the white 
dwarf in front of the M dwarf, an undetected infrared 
occultation observation and radial velocity observations 
of the M dwarf. 

2.2.4. Transit Light Curve 

We used a similar procedure to analyze the transit of 
the white dwarf in front of the M dwarf. Unlike the oc- 
cultation events, the transit events are not immediately 
apparent in the raw light curve due to the small geo- 
metric size of the white dwarf compared to the M dwarf. 
We de-trended and phase-folded the short-cadence light 
curve using the same procedure used in the previous sec- 
tion, except centered around the expected transit times 
assuming zero eccentricity: halfway between the occul- 
tation events. 

Agol (2003) showed that, for transits of compact ob- 
jects where the Einstein radius of the lens is significantly 
less than the radius of the source, the light curve obeys 
the same shape as a non-lensed planetary transit, but 
with a shallower or inverted transit depth. This is true 
even when limb- darkening of the source is included. This 
allows us to fit a conventional un-lensed transit light 
curve to the white dwarf transit data using an effec- 
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Fig. 4. — Long-cadence Kepler photometry of KOI-256. Each panel shows four days of relative flux, and panels are separated by several 
hundred days. All panels show a clear sinusoidal modulation, and the occultation events are seen as a brief drop in relative flux, each 
consisting of 2-3 long-cadence measurements. The transit events are not apparent in the light curve, and require phase-folding to detect 
(Section 2.2.4). The variability pattern shows the typical evolution in amplitude, phase, and overall shape, which is the signature of stellar 
activity. Compared to panel a, panel b shows variability with half the period and a smaller amplitude, and panel c shows variability at a 
different phase (note how the occultation is no longer close to time of minimum relative flux), and also a slightly smaller amplitude. We 
attribute this behavior to the combination of synchronous rotation of the M dwarf with the orbital period around the M dwarf-white dwarf 
center of mass, and an evolving spotted surface on the M dwarf which changes over 100 day timescales. 
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Fig. 5. — Signal residue vs. occultation period from the BLS anal- 
ysis for low- resolution and high- resolution near the peak(mser4). 
We use two methods to flnd the peak of the high-resolution power 
spectrum: a Gaussian-flt and spline-interpolation. Averaging the 
methods gives an occultation period of 1.3786548 ± 0.00001 days. 



live radius ratio parameter {R1/R2), which relates to the 
stars' physical radii and Einstein radius as such: 



(^\^ ^ ^WD - 2^Ein 
Vi?2 



p2 

-n-dM 



(1) 



where i?wD is the radius of the white dwarf, i?dM is the 
radius of the M dwarf and i?Ein is the Einstein radius. 
Assuming the source-to-lens distance is much smaller 
than the lens-to-observer distance, the Einstein radius 
takes the following form: 
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Fig. 6. — Phase-folded Kepler short-cadence light curve of KOI- 
256 during the occultation of the white dwarf as it passes behind 
the M dwarf. We have removed the sinusoidal modulation by fitting 
a polynomial to the out-of-eclipse data, normalized the resulting 
values, phased the data to match the period as measured using the 
BLS algorithm, and binned the data to 15 s windows. We fitted 
a model described in Section 2.2.3 to the un-binned light curve; 
however, we binned the data into 15 s windows in the figure for 
clarity. We show the residuals in the lower panel. 



where Mvvd is the mass of the white dwarf, a is the 
combined semi-major axis of the binary, G is Newton's 
gravitational constant and c is the speed of light. 

After subtracting the photospheric flux-contribution 
from the white dwarf using the occultation depth, we fit 
the transit model of Mandel & Agol (2002) to the data, 
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Fig. 7. — Phase-folded Kepler short-cadence light curve of KOI- 
256 during the transit of the white dwarf in front of the M dwarf. 
Solid line: best fitting model light curve. Dashed line: the model 
light curve expected without gravitational tensing. We applied the 
same detrending method used for the occultation events. We fitted 
a conventional planetary-transit model with limb-darkening from 
Mandel & Agol (2002) to the un-binned short-cadence light curve; 
however, we binned the data into 5-minute-wide windows in the 
figure for clarity. Since the M dwarf is significantly larger than the 
expected Einstein radius around the white dwarf during transit, 
the light curve follows the same shape as a planetary transit (Agol 
2003). 



using the expected quadratic limb-darkening coefficients 
of tfie M dwarf calculated by Claret & Bloemen (2011). 
We fixed the ingress/egress time and transit duration to 
match the values measured for the occultation, and fixed 
the mid-transit time to halfway in between the occulta- 
tion events. The only free parameters in the fit are the 
effective radius ratio {R1/R2) and a slope and offset for 
the out-of-eclipse light curve. We measure a combined ef- 
fective radius ratio of 0.0218 ± 0.0008. In Section 3.1, we 
use the effective radius ratio in combination with other 
measurements to determine the physical stellar masses 
and radii and the system orbital parameters. 

2.3. Infrared Occultation Light- Curve 

We observed an occultation event of KOI-256 with the 
Wide-Field Infrared Camera (WIRC) on the Palomar 
200-inch Hale Telescope (Wilson et al. 2003) on UT 2011 
August 19 from 04:50:53 UT to 08:12:29 UT. If the occul- 
tation were in fact a transiting planet, we would expect 
the same fractional depth in the infrared as in the Ke- 
pler bandpass. We used the _ff-band filter centered on 
1.635 /im and spanning ± 0.2 /im. We took 618 10 s ex- 
posures over the course of 202 minutes, corresponding to 
a duty cycle of 51%. We heavily defocused the telescope 
to spread out starlight and mitigate flat fielding errors, 
and "stared" on the target for the whole period. The 
airmass changed from 1.055 to 1.029 over the course of 
the observation. 

Reduction of the WIRC data was carried out using 
the routines developed in Zhao et al. (2012). We first 
subtracted the science images with corresponding aver- 
aged dark frames and a scaled "super-sky" image. The 
"super-sky" image was constructed with sky background 



frames taken right before and after the science frames. 
The data were then normalized with a master flat field. 
Thirteen reference stars within the flux range of 0.38 - 
1.1 times that of KOI-256 were selected to correct for the 
correlated common-mode systematics in the light curve. 
Aperture photometry was then performed on the tar- 
get and the reference stars, following the steps of Zhao 
et al. (2012). A final aperture with radius of 16.5 pix- 
els (4'.' 125) gave the smallest scatter in the data and was 
used for the photometry. A sky annulus with a 25 pixel 
inner radius and 35 pixel outer radius was used for back- 
ground estimation. Both the reduced out-of-occultation 
and in-occultation data closely follow the Gaussian noise 
expectations when binned (the scatter reduces as 1/VN, 
where N is the number of binning points), suggesting 
most of the correlated noise was corrected in the light 
curve. The final precision of the fitted depth of the light 
curve is 660 ppm, limited by the low brightness of the 
star and the systematics stemming from large centroid 
drifts (> 15 pixels) over the course of the observation. 

To determine the occultation depth, we fixed the tim- 
ing and duration to values determined in from the Ke- 
pler light curve and only allowed the depth to change 
in the fit. The UTC mid-exposure time of each images 
was converted to BJDtdb using the UTC2BJD routine 
by Eastman et al. (2010). The best-fit gives an occulta- 
tion depth of 0.00076 ± 0.00066, with a 3cr upper limit 
depth of 0.00275. The uncertainty was determined us- 
ing the largest value from the Levenberg-Marquardt fit, 
bootstrap, and residual permutation. Figure 8 shows the 
best-fit light curve model and the normalized data. 

The occultation depth in H-ha,nd is significantly 
smaller than the occultation depth in the Kepler band, 
consistent with the white dwarf occultation interpreta- 
tion. In Section 3, we calculate the stellar radii and 
white dwarf effective temperature (values listed in Ta- 
ble 3). For the radii and temperatures we measure, we 
expect an occultation depth of 0.001 in iJ-band, consis- 
tent with our upper limit, owing to the large flux contrast 
between the stars at this wavelength. 

2.4. Near-infrared Radial Velocities 

We observed KOI-256 with the TripleSpec Spectro- 
graph on the Palomar 200-inch Hale Telescope on UT 
2012 August 5, 6, 7, 8 and 11 for the purpose of measur- 
ing radial velocities of the M dwarf in KOT256. Triple- 
Spec is a near-infrared, long-slit spectrograph covering 
1.0 to 2.5 ^m simultaneously with a resolving power of 
2700 (Herter et al. 2008; Wilson et al. 2004). We ac- 
quired spectra of KOI-256 nine times over the course of 
five orbital periods to measure the object's radial veloc- 
ity variation. Each measurement consisted of four 80 s 
exposures, each taken at two positions along the Triple- 
Spec slit in order to subtract and remove contaminating 
sky emission lines. We reduced the spectra using the 
SpexTool program modified for TripleSpec at Palomar 
(Cushing et al. 2004, M. Cushing 2011, private communi- 
cation). For each observation of KOI-256, we acquired a 
spectrum of a nearby AOV star (HD 179095) to calibrate 
and remove telluric absorption lines using the xtellcor 
package (Vacca et al. 2003). Each resulting spectrum of 
KOI-256 had median signal-to-noise ratio of roughly 100 
per pixel in A'-band. 
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Fig. 8. — WIRC H-Hand light curve for a single occultation of 
the white dwarf by the M dwarf in KOI-256. We set a Scr upper 
limit of 0.00275 for the fractional occultation depth in H-band. 
The lack of a deep transit as seen in the Kepler data is consistent 
with the white dwarf interpretation of the system. In H-band, the 
expected flux ratio for the white dwarf is roughly 0.001, using the M 
dwarf and white dwarf temperatures and radius ratio determined 
in Section 3.2. 



The modified SpexTool program calculates a wave- 
length solution for reduced spectra using sky emission 
lines. We acquired deep (300 sec) exposures of the sky 
each night we observed KOI-256 to provide an accurate 
wavelength solution for each spectrum. We also observed 
G187 on UT 2012 August 6 using the same observing 
strategy, and achieved a similar S/N. G187 is a nearby 
(dis=10.4 pc, van Leeuwen 2007) M2.5 dwarf (Reid et al. 
1995) with an absolute radial velocity of 7.0 ± 2.0 km s~^ 
(Wilson 1953; Joy 1947). We measured an effective tem- 
perature of 3518 ± 50 K for G187 using the K-band in- 
dices of Rojas-Ayala et al. (2012). The effective tem- 
perature is similar to what was measured for KOI-256 
using the same indices {Tcs = 3450 ± 50 K, Muirhead 
et al. 2012). The similarity between the iiT-band spec- 
tra of the stars makes G187 a suitable target for use as 
a calibrated absolute radial velocity template spectrum 
for KOI-256. 

To measure radial velocities of the M dwarf component 
of KOI-256, we use a select region of iiT-band from 2.2 
to 2.4 /im, where M dwarfs have strong Na, Ca and CO 
features and the flux contribution from the white dwarf 
is negligible. First, we cross-correlated the spectrum of 
Gl 87 with each spectrum of KOI-256. The resulting 
radial velocity from the cross-correlation was used as a 
starting value for a Levenberg-Marquardt fit with three 
free parameters: two to adjust the continuum of the Gl 
87 spectrum to match KOI-256 and one to apply a radial 
velocity shift to the template. The Levenberg-Marquardt 
fit result was further refined and the posterior probability 
distribution explored using a Markov Chain Monte Carlo 
algorithm. The posterior probability distributions pro- 
vide more robust measurements of the true radial veloc- 
ity shifts, as the Levenberg-Marquardt fit is susceptible 
to local maxima in the distribution. The posterior prob- 
ability distribution also provides uncertainty estimates 



E 



o 
_o 
o 
> 

CO 

cc 



-150 

tn 

"§ 20r 

:g 0- 
s -20 r 




-100 



1244 1246 1248 1250 1252 
BJD - 2454900 

Fig. 9. — Radial velocities of the M dwarf in KOT256 measured 
with TripleSpec and using a spectrum of Gl 87 as a template and 
absolute RV reference, corrected for the motion of the observatory 
relative to the solar system barycenter. We fitted two Keplerian 
models to the data using the rvlin package (Wright & Howard 
2009): one with the eccentricity fixed to zero {solid line) and one 
with the eccentricity as a free parameter {dashed line, e = 0.056). 
We do not see a significant improvement in the fit; therefore we 
proceed assuming zero eccentricity. We fixed the orbital period 
and phase based on the occultation period and ephemeris as mea- 
sured using Kepler light curve (see Section 2.2.2, the occultation 
ephemeris times are marked with x's). The residuals are plotted for 
the zero-eccentricity model. The root-mean-square of the residuals 
is 8.4 kms~^, which is roughly twice the estimated per measure- 
ment uncertainty. 



for the measured radial velocities. The typical la uncer- 
tainty on the measurements is 4 km s~^. 

After correcting for the motion of the observatory rel- 
ative to the solar system barycenter and the archival ra- 
dial velocity measurement of Gl 87, we fitted a Keplerian 
model to the radial velocity measurements of KO-256 us- 
ing the rvlin package (Wright & Howard 2009). We 
fixed the orbital period and phase to that expected from 
the occultation events in the Kepler light curve. Figure 
9 shows our measured radial velocities and the Keplerian 
fit and Table 1 shows the best-fitting Keplerian parame- 
ters. The large amplitude and phase of the RV variations 
are consistent with a stellar-mass companion passing be- 
hind the M dwarf (occultation), rather than a planetary 
mass companion passing in front of the M dwarf (tran- 
sit, 180° phase difference). When combined with the 
trapezoidal shape and sharp ingress and egress of the oc- 
cultation light curve (Section 2.2.3), this provides strong 
evidence for the M dwarf-white dwarf interpretation of 
KOL256. 

2.5. High-resolution Optical Spectroscopy 

We observed the optical to far-red (3640-7820 A) spec- 
trum of KOI-256 using the High- Resolution Echelle Spec- 
trometer (HIRES Vogt et al. 1994). We used the stan- 
dard observing setup of the California Planet Survey 
(GPS; Howard et al. 2010; Johnson et al. 2010) includ- 
ing the red cross-disperser and the C2, 0'.'86 wide decker, 
with the iodine cell out of the light path. Our 500 s 
exposure resulted in an S/N of roughly 13 at 5500 A. 



TABLE 1 
TripleSpec Radial Velocities of the M Dwarf in KOI-256 

AND KePLERIAN FiT 



TABLE 2 
Equivalent Width Measurements for KOI-256 



BJD - 2454900 


RV (km s-i) 


1244.6636 


96.1 ± 3.6 


1245.6961 


6.8 ± 5.4 


1245.8384 


69.0 ± 4.5 


1245.9457 


85.1 ± 3.9 


1246.7308 


-116.1 ± 3.6 


1246.9384 


-80.1 ± 3.9 


1247.7754 


-6.4 ± 3.3 


1250.6623 


-79.6 ± 3.7 


1250.8408 


-109.0 ± 3.7 


Parameter 


Value 


Per (days, fixed) 


1.3786548 ^ 


Tp (BJD, fixed) 


2455374.324838^ 


e (fixed) 


0.0 


uj (fixed) 


90.°0 


K (km s-i) 


105.9 ± 1.6 


7 (km s"-*-) 


-7.3 ± 2.0'' 


dv/dt (fixed) 


0.0 



Tixed to the oceultation period and ephemeris measured from 
the Kepler Hght eurve (see Seetion 2.2). 

''Measured relative to Gl 87, assuming literature value of 7.0 it 
2.0 km s-i (Wilson 1953; Joy 1947). 



2.5.1. Vsini of the M dwarf 

Figure 10 shows the HIRES speetra of KOI-256 and a 
nearby M dwarf monitored by CPS, Gl 179, in molecular 
bands corresponding to select Palomar Michigan State 
University (PMSU) indices including Ha (Reid et al. 
1995). Gl 179 is an M3 Dwarf with a measured effective 
temperature Toff =3424 ± 16 K and overall metallicity 
[M/H] = -f-0.17 ± 0.12 dex as measured by Rojas-Ayala 
et al. (2012). Muirhead et al. (2012) used the same K- 
band spectroscopic techniques to measure similar values 
for KOI-256 (T^s = 3450 ± 50K, M/H = -hO.31 ± 0.10). 
Browning et al. (2010) measured the rotational broad- 
ening of Gl 179 to be Vsini < 2.5kms~^ using archival 
CPS data. With such a low Vsini and similar Tgg , 
Gl 179 is ideal for use as a template for measuring the 
V sini of KOI-256. We combined all of the spectra taken 
of Gl 179 over the course of the CPS survey using regions 
that are not contaminated by molecular iodine absorp- 
tion. The result is a spectrum with an S/N of roughly 
2500 per pixel. The large radial velocity change of the M 
dwarf as it orbits the white dwarf results in a smearing 
of the HIRES spectrum during the exposure. However, 
when added in quadrature to the rotational broadening, 
the smearing effect contributes less than our estimated 
Vsini measurement errors. 

To measure the Vsini of KOI-256, we following the 
approach of Browning et al. (2010). For each spectral 
region in Figure 10 (excluding Ha) we cross-correlated 
the high S/N spectrum of Gl 179 with a rotationally 
broadened version of itself for a grid of V sin i values 
spanning 1 to 30 km s^^. We then fit a Gaussian pro- 
file to the resulting cross-correlation function (XCFs), 
and recorded the fuU-width-at-half maximum: the XCF 
FWHM. With the resulting calibration between XCF 
FWHM and Vsini , we cross-correlated the unbroad- 



Line 


EW (A) 


Line Region (A) 


Continuum (A) 


Ha 


3.40 ± 0.06 


6561-6565 


6557-6560, 6566-6569 


H7 


3.87 ± 0.26 


4339-4342 


4335-4338, 4344-4347 


m 


3.67 ± 0.42 


4100-4104 


4096-4099, 4105-4108 


He 


3.97 ± 0.21 


3969.5-3971 


3971-3975 


HC 


6.02 ± 0.48 


3887-3890 


3883-3886, 3892-3895 


CallH 


16.28 ± 3.33 


3968-3969.5 


3971-3975 


CallK 


13.56 ± 0.46 


3932-3934 


3928-3930, 3936-3939 



ened Gl 179 spectrum with the KOI-256 spectrum, and 
recorded the XCF FWHM for each spectral region. We 
then interpolated the XCF FWHM onto the calibration 
XCF FWHM versus V sin i relation. 

The procedure requires an assumption about the lin- 
ear limb-darkening coefficient for the broadening kernel 
used on the Gl 179 template spectrum. We use the 
value 0.6715 calculated by Claret & Bloemen (2011), cor- 
responding to the limb-darkening expected from an M 
dwarf with T^s =3400 K and log{g) = 5.0. We take 
the weighted mean of the measurements for the bands as 
our best estimate (19.67 km s^^) with an uncertainty of 
0.52 km s~^. Figure 11 illustrates the consistency of the 
measured V sin i for each of the different bands and the 
weighted mean. 

2.5.2. Emission Line Strengths 

We note that KOI-256 shows strong Balmer emission 
lines and Call H & K. Table 2 lists equivalent widths of 
selected emission lines as measured in the HIRES spec- 
trum, and Figure 10 shows the strong Ha emission rel- 
ative to Gl 179. We also report the line and continuum 
spectral regions used to calculate the equivalent widths. 
We calculated the uncertainty in the equivalent widths 
by injecting noise into the spectra and recalculating. The 
emission likely comes from chromospheric activity in the 
M dwarf related to rapid rotation. This is consistent 
with the rotation-activity correlation for PCEB M dwarfs 
measured by Morgan et al. (2012) using spectra from the 
Sloan Digital Sky Survey. 

2.6. Low-resolution Optical Spectroscopy 

We observed KOI-256 with the Echellctte Spectro- 
graph and Imager (ESI) on the Keck II Telescope (Sheinis 
et al. 2002) on UT 2012 October 9, for the purpose of de- 
tecting the contribution to the spectrum from the white 
dwarf companion at blue wavelengths. We operated ESI 
in "low-D" mode, using a I'.'O slit to obtain a single-order 
spectrum from 3900 Ato 11000 A. The resolving power 
of ESI in low-D mode varies from R=6000 at 3900 Ato 
R=1000 at 11000 A. We co-added 5 separate 300 s expo- 
sures, so as to acquire enough flux at the blue end of the 
detector without saturating the red end, which contains 
significant flux from the M dwarf component. 

We also observed the spectro-photomctric standard 
BD+28 421 using identical settings but with a single ex- 
posure of 10 s of integration.^^ After subtracting a bias 
frame from the raw images and dividing by a high signal- 
to-noise twilight flat, we extracted both stars' spectra us- 

^^ http;//www. eso.org/sci/observing/tools/standards/spectra.html 
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Fig. 10. — HIRES spectra of KOI-256 (black), Gl 179 (blue) and a rotationally broadened version of Gl 179 (broadened to Vsini = 19.6 
km s"'^, red) for select PMSU molecular bands (Reid et al. 1995) including Ha. The spectra are continuum normalized and manually 
shifted in wavelength to overlap within each spectral region. Gl 179 has similar Teg and [M/H] to the M dwarf component of KOI-256, as 
measured by the Rojas-Ayala et al. (2012) methods. The KOI-256 spectra are dominated by the M dwarf component, which constitutes 
nearly 98% of the flux in these regions based on the occultation depth. The M dwarf in KOI-256 is significantly rotationally broadened 
compared to to Gl 179, a known slow rotator (Browning et al. 2010). We measure the V sini of the M dwarf to be 19.6 ± 1.0 km s~^ 
(see Figure 11). We note the strong Ho emission in the M dwarf, for which we measure an equivalent width of 3.32 ± 0.05 A. The Ha 
emission likely comes from chromospheric activity in the rapidly rotating M dwarf, consistent with the rotation-activity correlation in M 
dwarfs (Morgan et al. 2012). 
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Fig. 11. — Plot showing the measured V sini of the M dwarf in 
KOI-256 using different feature-rich regions of the HIRES spectrum 
corresponding to the bands designated TiO 1, TiO 3, TiO 4, CaH 
2 and CaH 3 by Reid et al. (1995), as well as the mean from all 
regions. We follow the procedure of Browning et al. (2010) to mea- 
sure Vsinj , using Gl 179 as a slowly-rotating template spectrum 
(see Section 2.5). We take the weighted mean of the measurements 
as our best estimate (19.67 km s~^) with an uncertainty of 0.52 
km s~^. 



ing the ATV image display software (Barth 2001). We fit- 
ted a custom wavelength solution to the extracted spec- 
tra based on the positions of selected, bright HgXe emis- 
sion lines in a calibration exposure. Figure 14 shows the 
flux-calibrated ESI spectrum. The optical spectrum is 
consistent with an M3 dwarf star, and we include a tem- 
plate spectrum from Pickles (1998) in the figure. We do 
not, however, detect any flux from the white dwarf in 
the blue part of the spectrum. In Section 3.2, we convert 
the M3 template spectrum to a surface brightness using 
an empirical radius measurement of another M3 star, Gl 
581, for the purpose of determining the white dwarf ef- 
fective temperature. We show that the white dwarf is 
relatively cool and therefore should not appear in the 
ESI spectrum. 

2.7. GALEX NUV Photometry 

KOI-256 was observed with the near-ultraviolet (NUV) 
channel on the Galaxy Evolution Explorer (GALEX, 
Martin et al. 2003) spacecraft on UT 2012 August 7, 
8, 13, 16, 25, 28, 29 and 31, and 2012 September 02, 05, 
08 and 09 as part of an independent program to mosaic 
the Kepler field (J. Lloyd 2012, private communication). 
The NUV band in GALEX spans roughly 1000 Awith 
an effective wavelength of 2271 A. Figure 12 plots the 
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resulting flux in /iJy versus time and orbital phase. 

The observation at BJD 2456155.9617973 represents a 
6(7 increase in the NUV flux from KOI-256 and may be 
a flare event. The slow increase in NUV flux from BJD 
2456165 to 2456180 may be due to general increase in 
chromospheric activity. With the release of quarter 14 
data from the Kepler Mission, we can correlate the NUV 
and visible light curves. One measurement lies within 
the occultation and one within the transit, however nei- 
ther is significantly different than the out-of-eclipse ob- 
servations. Given the strong Ha emission from the M 
dwarf, and the low white dwarf temperature we measure 
in Section 3.2, we ascribe the NUV flux to chromospheric 
activity in the M dwarf and not the white dwarf compo- 
nent. 

3. SYSTEM PARAMETERS 

In this section, we combine the measurements to em- 
pirically determine the M dwarf and white dwarf param- 
eters, which are listed in Table 3. 

3.1. Masses, Radii and Orbital Parameters 

To determine the stellar masses, radii and system or- 
bital parameters, we combine six measurements: the or- 
bital period (P), the semi-amplitude radial velocity of 
the M dwarf (Kdu), the occultation duration (idur), the 
ingress/egress duration (iing), the effective radius ratio 
parameter fitted to the transit light curve {R1/R2) and 
the rotational broadening of the M dwarf {Vsini )-with 
six unknowns-the M dwarf mass {M^m), the M dwarf 
radius (RdM), the white dwarf mass (Mwd)) the white 
dwarf radius (i?wD), the orbital inclination (i) and the 
system semi-major axis (a)-in six equations. The first 
equation is Equation 1 in Section 2.2.4, relating the ef- 
fective radius ratio of the transit light curve to the stellar 
radii, white dwarf mass, and semi-major axis. The sec- 
ond equation is Kepler's Third Law, which relates the 
measured orbital period to the masses and semi-major 
axis: 



P = 27r. 



G (MdM + A^wd) 



(3) 



The final four equations are described below. 

Assuming spin-orbit synchronization, as evidenced by 
the out-of-eclipse variations in the long-cadence light 
curve (see Section 4 and Figure 4) , the rotational period 
and spin inclination of the M dwarf match the orbital 
period and orbital inclination of the binary. In this case, 
the measured V sin i of the M dwarf relates to the or- 
bital period, orbital inclination and M dwarf radius by 
the following equation: 



V sin i ~ 



2tt RdM sin i 
P 



(4) 



The measured semi-amplitude radial velocity of the M 
dwarf relates to the stellar masses and orbital period by 
the mass function. Assuming zero eccentricity: 



K. 



dM 



27r G \ 3 MwD sin i 



(5) 



P ' (MwD + MdM)3 

The measured occultation duration (mid-ingress to 



mid-egress) relates to the system parameters by the fol- 
lowing equation: 



tdx. 



2 i?dM sin i Vl — fo^ 
Ka^ (1 + Mam/M^-d) 



(6) 



where b = a cosi/Rdu is the impact parameter of the 
white dwarf relative to the center of the projected M 
dwarf on the sky, in units of the M dwarf radius. And 
finally, the measured ingress/egress time (ting) follows 
the following equation, treating the limb of the projected 
image of the M dwarf as a straight edge, which the white 
dwarf passes across at non-normal incidence: 



tu 



2 RdM sinz 



KdM (1 + MdM/Mwo) VT^^ 



(7) 



We solve for the six unknowns using Broyden's 
Method, a secant method for solving nonlinear simulta- 
neous equations (Broyden 1965). To estimate the uncer- 
tainties in the solved parameters due to errors in the mea- 
surements, we take a Monte Carlo approach. We created 
1000 copies of the measured values and added normally- 
distributed noise corresponding to the measurement er- 
rors to each copy. For each noise-added copy of the mea- 
sured parameters, we then solved for the unknowns, and 
take the standard deviations in the resulting 1000 iter- 
ations as the Ict uncertainties. We report the resulting 
stellar masses, radii, orbital parameters and correspond- 
ing uncertainties in Table 3. 

3.2. White Dwarf Effective Temperature 

The white dwarf is too cool to detect in the blue wave- 
lengths of the Keck ESI spectrum (Figure 14), and the 
GALEX measurements are likely contaminated by UV 
emission from the active M dwarf chromosphere. We 
choose to measure the white dwarf temperature using the 
occultation depth in the Kepler band. The fractional oc- 
cultation depth (Sc) equals the ratio of the flux from the 
white dwarf to the total flux from the system within the 
Kepler bandpass: 



5c = 



f 



Tkop.A SwD,x Rwn ^ d\ 



J (Tkcp.A S'wD,A -R-WD + ^Kcp,A SdM.A ^dlvl) ^ '^^ 



(8) 



Where Tksp.a is the transmission of the Kepler pho- 
tometer versus wavelength, S'wd.a is the surface bright- 
ness of the white dwarf per unit wavelength, and S'dM,A 
is the surface brightness of the M dwarf per unit wave- 
length. The surface brightnesses of the two bodies de- 
pend on their respective temperatures and opacities over 
a wide range of wavelengths. With a model for the sur- 
face brightness of the M dwarf that is consistent with the 
measured ESI spectrum, and a temperature-dependent 
model for the surface brightness of the white dwarf, we 
can vary white dwarf temperature to match the occulta- 
tion depth in the Kepler bandpass. 

For the M dwarf, synthetic surface brightnesses are a 
matter of signiflcant debate. Model atmospheres have 
historically been unable to match spectroscopic observa- 
tions due to the complicated molecular opacities in their 
atmospheres. The gold standard has been the PHOENIX 
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Fig. 12. — GALEX NUV flux from KOI-256. Left: flux in /ijy vs. barycentric Julian date (BJD) of the observation. Right: same, 
but phase folded to the occultation period with the occultation ingress and egress ephemera represented by dashed lines and the transit 
by dotted lines. Given the strong Ho emission in the M dwarf and low white dwarf temperature (see Section 3.2), we ascribe the NUV 
emission to the M dwarf chromosphere and not the white dwarf photosphere. 



model atmosphere code (e.g. Hauschildt et al. 1999), 
which reproduce infrared spectra of M dwarfs relatively 
well (Rojas-Ayala et al. 2010, 2012). However, at optical 
wavelengths the spectra show discrepancies from obser- 
vations (e.g. Allard et al. 2012). 

Instead of using a model surface brightness for the 
M dwarf, we chose to use template spectra from Pick- 
les (1998) that have been flux-calibrated to match M 
dwarfs with accurate angular diameter measurements us- 
ing optical long-baseline interferometry. Boyajian et al. 
(2012) measured angular diameters for GJ 411 (M3V) 
and GJ 699 (M4V) and von Braun et al. (2011) mea- 
sured the angular diameter of Gl 581 (M3V), both using 
the CHARA Array on Mt. Wilson and correcting for 
effects from limb-darkening. In their respective papers, 
they flux-calibrated M2, M3 and M4 template spectra 
measured by Pickles (1998) to match the visible and in- 
frared photometry of these stars. By dividing the flux 
calibrated template spectra by the measured angular di- 
ameters measured for the star, we create semi-empirical 
surface brightness templates for the M dwarf in KOI-256. 
We use these three stars, as they have similar spectral 
type to the M dwarf in KOI-256. 

For the white dwarf, we used the Planck function to 
model the surface brightness. To determine the white 
dwarf temperature, we computed the expected occulta- 
tion depth as a function of the Planck function temper- 
ature in Equation 8. The absolute value of the differ- 
ence between the calculated and measured occultation 
depth was used as a goodness-of-fit statistic. The best 
fitting temperatures were found to be 6100 K, 7100 K 
and 7400 K using the M2, M3 and M4 surface brightness 
templates for the M dwarfs respectively (see Figure 13). 
We chose the value 7100 K as the best estimate because 
the spectral type of the M dwarf in KOI-256 is closest 
to M3, and take the standard deviation in the three es- 
timates as the uncertainty (700 K). Figure 14 plots the 
measured spectrum of KOI-256 with the model spectra, 
Kepler bandpass and GALEX flux measurements. 



4. DISCUSSION 
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Fig. 13. — Absolute value of the observed minus calculated Ke- 
pler fractional occultation depth versus assumed white dwarf effec- 
tive temperature. We use the Planck function to model the white 
dwarf surface brightness, and the flux calibrated spectra and an- 
gular diameters of Gl 411 (M2, dot-dashed line), Gl 581 (M3, solid 
line) and Gl 699 (M4, dashed line) from von Braun et al. (2011) and 
Boyajian et al. (2012) to model the M dwarf surface brightness. We 
solve for the white dwarf effective temperature numerically, using 
the observed minus calculated Kepler fractional occultation depth 
as a goodness-of-fit parameter. 



KOI-256 is one of several eclipsing systems containing 
a white dwarf discovered in the Kepler field. Four of 
them, KOI-74 and KOI-81 (Rowe et al. 2010; van Kerk- 
wijk et al. 2010), KIC 10657664 (Carter et al. 2011), and 
KOI-1224 (Breton et al. 2012) are composed of a low- 
mass, hot white dwarf, with Mwd ~ 0.2Mq and Tgg > 
10,000 K, and an early-type main sequence star. Two 
other systems have a hot WD whose mass is close to the 
canonical mass of 0.6 Mq. Those are KPD 1946-f4340 
(Bloemen et al. 2011), where the WD is orbiting a sub- 
dwarf B star, and KIC 10544976 (Almenara et al. 2012), 
with an M dwarf primary. The latter two have an orbital 
period shorter than 10 hrs, compared to the 33 hr period 
of KOI-256. Therefore, KOI-256 is the first cold WD dis- 
covered in the Kepler field orbiting a low-mass star at a 
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Fig. 14. — Observed and modeled spectra of KOI-256. The Keck-ESI spectrum {solid black line) is especially noisy from 4000 to 5000 A. 
The model spectra are computed using the surface brightness of an M3 dwarf, calibrated using the angular diameter measurement of Gl 
581 (von Braun et al. 2011) and a Plank function with a temperature of 7100 K for the white dwarf. We then scaled the surface brightness 
to the expected flux using the stellar radii measurements described in Section 3.1 and KOI-256's r-band magnitude in the Kepler Input 
Catalog (KIC, Batalha et al. 2010). The GALEX flux measurements are inconsistent with the photosphere of cither the M dwarf or the 
white dwarf, and are likely enhanced by the chromosphere of the M dwarf. 



relatively long period. KOI-256 is also a model example 
of gravitational lensing effects during transit. With em- 
pirically measm'cd masses and radii, KOI-256 serves as a 
benchmark object for studying the effects of PCEB evo- 
lution on M dwarf abundance and activity as a function 
of mass and radius. 



4.1. The Evolution of the System 

With its short orbital period and white dwarf sec- 
ondary, KOI-256 is unquestionably a PCEB and will 
eventually evolve into a cataclysmic variable. However, 
KOI-256 has a relatively long orbital period compared 
to most known "pre-CVs" (Schreiber & Gansicke 2003), 
and it contains a relatively cool white dwarf. These char- 
acteristics indicate KOI-256 is relatively old. The time 
elapsed since the end of the common envelope phase can 
be estimated from the white dwarf cooling time. Using 
the cooling tracks of Renedo et al. (2010), we estimate 
a cooling time icooi ~ 2 Gyr for the white dwarf in KOI- 
256. We also use the initial to final mass relationships 
of Renedo et al. (2010) to estimate that the white dwarf 
in KOI-256 originated from a progenitor of initial mass 
1.5M0 < Mi < 3Mq, although the initial to final mass 
relationship is less certain for white dwarfs in PCEB bi- 
naries (Rebassa-Mansergas et al. 2011). 

The orbits of pre-CVs decay primarily due to angular 
momentum loss from the dwarf star's magnetized wind. 



We use the angular momentum loss prescription of Sills 
et al. (2000) and Andronov (2003), in which the system's 
angular momentum, J, changes as 



dJ 



^ = -K 



M 



1/2 



R 

Rq 



1/2 



W^ U) 



(9) 



where K = 2.7 x lO'^^g cm^ s is a calibrated constant, 
uj is the angular spin frequency of the dwarf star, and 
ujc is the angular spin frequency at which the angular 
momentum loss saturates (we use Wc = Twq). Assuming 
efficient spin-orbit coupling and ignoring the small angu- 
lar momentum contained in the stars, the initial orbital 
frequency of the system was 



n. 



AAt. 



cool 



a 



-4/3 



-1-3/4 



(10) 



where /i is the reduced mass of the system, Mt is its 
total mass, A = K[RMq/{MRq)Y/^oj'^, and Qq is the 
current orbital angular frequency. We use equation (10) 
to estimate that KOI-256 emerged from the common en- 
velope phase with an orbital period of Pi « 1.7 days. 
The elapsed time until Roche lobe overflow is 



tRL 



KGM,)y^ 

AA 



(n- 



4/3 






(11) 



where ^Irl is the angular orbital frequency at Roche lobe 



13 



TABLE 3 

Measured and Calculated Parameters for KOI-256 



Parameter 


Value 


Notes and Relevant Section 


Spectroscopic Parameters 










M dwarf semi-amplitude RV (i^dM) 


106.5 ± 1.8 km s-i 






2.4 


Orbital eccentricity 


< 0.05 


Assumed to be 0.0, 2.4 


M dwarf V sin i 


19.79 ±0.52 km s"! 






2.5 


Occultation Parameters 










Orbital period [days] 


1.3786548 ± 0.000010 






2.2.2 


Mid-occultation cphemeris (BJD) 


2455373.635498 ± 0.000036 






2.2.2 


Occultation duration (tdur) 


60.051 ± 0.037 minutes 






2.2.3 


Ingress/egress duration (ting) 


1.604 ± 0.067 minutes 






2.2.3 


Kepler fractional depth (5c ) 


0.0238959 ± 0.0000047 






2.2.3 


H-hand fractional depth 


< 0.00275 


3o- 


upper 


limit, 2.3 



Transit Parameters 

Effective radius ratio {R1/R2) 

M dwarf Parameters 

Mass (MdM) 
Radius (-Rdlvl) 
Effective temperature 
Metallicity [M/H] 
Kepler magnitude (Kp) 
Kr,-H 



0.0209 ± 0.0011 



0.51 ± 0.15 Mq 

0.540 ± 0.014 Rq 

3450 ±bO K 

+0.31 ± 0.10 

15.40 

3.37 



2.2.4 



3.1 
3.1 

from Muirhead et al. (2012) 

from Muirhead et al. (2012) 

KIC and occultation depth 

KIC 



White D^varf Parameters 

Mass (Mwd) 

Radius (-Rwo) 

Einstein radius (i?Ein) 

WDTeff 

Kepler magnitude {Kp) 

Kp-H 

Orbital Parameters 

Combined semi-major axis (a) 
Orbital and M-Dwarf-Spin Inclination 

a/RdM 

Impact parameter (b) 



0.592 ± 0.084 Mq 

0.01345 ± 0.00091 Rq 

0.00473 ± 0.00055 Rq 

7100 ± 800 K 

19.45 

< 1.05 



0.0250 ± 0.0018 AU 

89.°01 ± 0.65 

9.97 ± 0.64 

0.17 ± 0.10 



3.1 
3.1 
3.1 

3.2 

KIC and occultation depth 
H-hanA occultation upper limit 



3.1 
3.1 
3.1 
3.1 



Assumptions 

The M dwarf is synchronously rotating with orbit. 

The M dwarf spin axis is aligned with binary orbital axis. 

The M dwarf obeys Claret & Bloemen (2011) limb-darkening parameters. 



Evidence in 2.2.1 
Evidence in 2.2.1 



overflow. We estimate mass transfer will commence at 
an orbital period of Phl — 280 minutes in KOI-256, and 
that this will occur after tji^ ~ 5.3 Gyr of orbital decay. 
After the initiation of mass transfer, KOI-256 will likely 
evolve as a fairly typical CV system. 

We can also calculate whether the orbital decay in 
KOI-256 is likely to be observed via eclipse timing. Over 
an observational baseline T, the eclipse times change by 
an amount 

M=-^. (12) 

In KOI-256, we estimate the eclipse times will only 
change by At « 1 s over the next twenty years, so it 
is unlikely that the orbital decay will be detected in the 
near future. 

4.2. On the High Metallicity of the M Dwarf 

Muirhead et al. (2012) measured an overall metallic- 
ity of [M/H] = +0.31 ± 0.10 using the K-hanA near- 



infrared spectroscopic indices of Rojas-Ayalaet al. (2010, 
2012). The indices are calibrated on widely separated 
main-sequence binary stars, and the calibration may not 
be applicable to the M dwarf in KOI-256 due effects 
from common-envelope evolution. Nevertheless, the high 
metallicity compared to the other M dwarf KOIs ana- 
lyzed in that paper (mean [M/H] = -0.11, RMS [M/H] 
= 0.15) is interesting considering that nearby M dwarfs 
with white dwarf companions also show high metallicities 
using the K-hand calibration. Rojas-Ayala et al. (2012) 
found that G 203-47 and GL 169.1 A, both nearby M 
dwarfs with white dwarf companions, had metallicities of 
[M/H] = 0.24 ± 0.12 and 0.26 ± 0.12, respectively. They 
speculate that the high metallicities could be due to en- 
hancement by nucleosynthesis products from the white- 
dwarf progenitor during the common-envelope phase. If 
true, the high metallicity of M dwarfs in PCEBs maybe 
a concern for metallicity-biased exoplanet surveys of M 
dwarfs since their high metallicity is not primordial and 
therefore does not necessarily indicate an increased like- 
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lihood of planet occurrence (e.g. Apps et al. 2010). 

4.3. Future Observations 

Future investigation of KOI-256 will enable ever more 
precise measurements the parameters of the system, in- 
cluding the mass, radius, and temperature of both ob- 
jects, which are especially interesting in the case of 
the white dwarf. These investigations will benefit from 
follow-up observations beyond future Kepler data, in- 
cluding UV spectroscopy and multi-band high-speed oc- 
cultation photometry. The latter will constrain the radii 
and temperature ratio, and may be able to detect the 
asymmetry in the occultation ingress and egress light 
curves due to the Photometric Rossiter-McLaughlin ef- 
fect (Shporer et al. 2012; Groot 2012). Modehng this 
effect will enable measuring the white dwarf obliquity 
and rotational velocity. 
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